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ABSTRACT 

Context. The dynamo processes in cool active stars generate complex magnetic fields responsible for prominent surface stellar 
activity and variability at diff'erent time scales. For a small number of cool stars magnetic field topologies were reconstructed from 
the time series of spectropolarimetric observations using the Zeeman Doppler imaging (ZDI) method, often yielding surprising and 
controversial results. 

Aims. In this study we follow a long-term evolution of the magnetic field topology of the RS CVn binary star II Peg using a more 
self-consistent and physically more meaningful modelling approach compared to previous ZDI studies. 

Methods. We collected high-resolution circular polarisation observations of II Peg using the SOFIN spectropolarimeter at the Nordic 
Optical Telescope. These data cover 12 epochs spread over 7 years, comprising one of the most comprehensive spectropolarimetric 
data sets acquired for a cool active star. A multi-line diagnostic technique in combination with a new ZDI code is applied to interpret 
these observations. 

Results. We have succeeded in detecting clear magnetic field signatures in average Stokes V profiles for all 12 data sets. These profiles 
typically have complex shapes and amplitudes of ~10~^ of the unpolarised continuum, corresponding to mean longitudinal fields of 
50-100 G. Magnetic inversions using these data reveals evolving magnetic fields with typical local strengths of 0.5-1.0 kG and 
complex topologies. Despite using a self-consistent magnetic and temperature mapping technique, we do not find a clear correlation 
between magnetic and temperature features in the ZDI maps. Neither do we confirm the presence of persistent azimuthal field rings 
found in other RS CVn stars. Reconstruction of the magnetic field topology of II Peg reveals significant evolution of both the surface 
magnetic field structure and the extended magnetospheric field geometry on the time scale covered by our observations. From 2004 
to 2010 the total field energy drastically declined and the field became less axisymmetric. This also coincided with the transition from 
predominantly poloidal to mainly toroidal field topology. 

Conclusions. A qualitative comparison of the ZDI maps of II Peg with the prediction of dynamo theory suggests that the magnetic 
field in this star is produced mainly by the turbulent dynamo rather than the solar aQ dynamo. Our results do not show a clear 
active longitude system, nor is there an evidence of the presence of an azimuthal dynamo wave. 
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• 1. Introduction 
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Direct, spatially -resolved observations of the solar disk demon- 
strate the ubiquitous presence of magnetic fields and their key 
role in driving time-dependent, energetic processes. It is be- 
lieved that magnetic fields are generated in the solar interior, 
due to a complex interplay between the non-uniformities in the 
internal rotation profile, large-sc ale flows, and vigo rous turbu- 
lence due to convective motions ( |Ossendrijver|2QQ3| ). Although 
the rotation profile, with two prominent regions of high radial 
shear near the bottom and top of the solar convection zone, 
is known rather reliably from helioseismic observations (e.g. 
[Thompson et al. 2003), measurements of the potentially very 
important meridional flow pattern, e specially deeper in the con - 
vection zone, are highly uncertain ( [Gough & Hindman||2010|). 
A variety of models of the solar dyna mo exists (e.g. [Dikpati & 



2012| ), but none of them is capable of fully reproducing the ob- 
served characteristics of the solar magnetic cycle. Direct numer- 
ical simulations aiming at including all the relevant physics in a 
single model, with the expectation being that all the observable 
features from the flow patterns up to the magnetic cycle would 
emerge self-consistently, have so far been only marginally suc- 
cessful ( [Ghizaru et aT|[2QTQ[ [Racine et al.]|25TT| jKapyla et al^ 
2012a| ). Therefore, it is fair to say that the solar cycle remains 



Charbonneau[1999[[Kapyla et al.[2006t[Kitchatinov & Qlemskoy 



^ Based on observations obtained with the Nordic Optical Telescope, 
operated on the island of La Palma jointly by Denmark, Finland, Ice- 
land, Norway, and Sweden, at the Spanish Observatorio del Roque de 
los Muchachos of the Instituto Astrofisica de Canarias. 



a theoretical challenge, with some key information on the large- 
scale flow patterns yet to be supplied by observations. 

Compared to the Sun, many cool stars exhibit enhanced ac- 
tivity levels, suggesting the presence of stronger surface mag- 
netic fields produced by more eflftcient dynamos. Observations 
of such active stars oflfer possibilities to constrain theoretical 
models by examining the eflftciency of dynamo mechanisms as a 
function of stellar mass, age, and rotation rate. However, direct 
magnetic field analysis of cool active stars is impeded by the 
weakness of polarisation signatures in the spatially-unresolved 
stellar spectra. Furthermore, dynamos in cool stars usually give 
rise to topologically complex magnetic fields which cannot be 
meaningfully characterised using observational techniques sen- 
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sitive only to the global m agnetic field component ([ Borra & 
Landstreetl 19801 |Vogt| 198Q| ). 



A major breakthrough in studying the magnetism of cool 
stars was achieved by combining multi-fine polarisation diag- 
nostic techniques ( [Donati et al.| 1997 ) applied to high-resolution 
spectropolarimetric observations with the inversion procedure of 
Zeeman Doppler imaging (ZDI, |Brown et al.|| 19911 ponati & 



|Brown 1997). This approach, capable of providing spatially- 
resolved information about stellar magnetic field topologies, was 
applied to a number of late-type stars (e.g. Dona ti et"aL]|2003 



Petit et al.|2QQ 4), but only three objects - AB Dor, LQHya, and 



HR 1099 - were systematically followed with multiple magnetic 
maps on the time scales of five to ten years. 

The results of ZDI mapping are often spectacular in their 
ability to resolve fine details of stellar magnetic field geome- 
tries, but also controversial in several aspects. First, analysis of 
stellar spectra is limited to only circular polarisation because, be- 
ing roughly ten times weaker, the linear polarisation signatures 
are not readily detectable at the moderate signal-to-noise (S/N) 
ratios used for the circular polarisation monitoring ( |Kochukhov| 
|et al.|[20TT] ). Supplying inversion algorithms with such incom- 
plete Stokes vector data is known to lead to spurious features 
in the reconstructed magnetic maps ( [Donati & Brown||1997| 
|Kochukhov & Piskunov^ 2002 ). Second, interpretation of the 
mean-line shapes constructed by averaging thousands of lines is 
inevitably inferior to modelling individual spectral lines because 
of the d iflSculty in choo sing appropriate mean-line parameters 
(Kochuk hov et al.||2010| ). Finally, and most importantly, nearly 
all magnetic field maps reconstructed from the Stokes V observa- 
tions of cool active stars were obtained separately and inconsis- 
tently from the mapping of brightness distributions from Stokes 
/, raising questions about the validity of magnetic maps and pre- 
venting direct analysis of the spatial relation between magnetic 
and temperature inhomogeneities. 

These considerations suggest that substantial progress in un- 
derstanding cool-star magnetism through ZDI requires continu- 
ing research in several directions. On the one hand, the sam- 
ple of active stars studied with multiple-epoch magnetic images 
has to be expanded to justify the far-reaching conclusions previ- 
ously inferred from the analysis of a few objects. On the other 
hand, ZDI studies must re-examine key methodological limita- 
tions of this technique and strive to employ physically realistic 
modelling approaches whenever possible. This paper addresses 
both of these aspects in a detailed ZDI analysis of the RS CVn 
binary II Peg. 

The source of the most powerful stellar flares ever observed 
( Osten et al. 2007 ) and the brightest X-ray object within 50 pc 
fMakarov 2003 ), II Peg (HD 224085, HIP 117915) is one of the 
most prominent active cool stars in the solar neighbourhood. 
This star is an RS CVn-type, single-line spectroscopic binary 
(SB 1) with an orbital period of ^ 6.72 days, consisting of a K2IV 
primary and a low-mass (M0-M3V) secondary. The primary 
star exhibits copious manifestations of the magnetically-driven 
surface activity, including a strong non-thermal emission in the 
UV and optical chromospheric lines and in the X-ray and ra- 
dio wavelength regions. It also shows powerful flares as well as 
regular pho tometric and spectroscopi c variations due to evolving 
cool spots. [Berdyugina et al.| ( [1998b) studied the orbital motion 
of the massive component and provided a comprehensive sum- 
mary of the physical properties of II Peg. Many studies exam- 
ined photometric variations of this star (e.g. Siwak et al. 2010; 



the role of active longitudes ( [Berdyugina et al.||l999t [Rodono 
[et al.|2000| . 

Other studies targeted II Peg with high-resolution spectro- 
scopic observations with the goal to constrain the spot tempera- 
tures using TiO absorption bands ( [O'Neal et al.| 19981 [Berdyug-] 
[ina et al.|1998b |) and to analyse configurations of the surface tem- 
perature inhomogeneities with the Dopp l er im aging technique 
([Berdyugina et al.,1998al[Gu et al.|2003t[Lindborg et al.|2011 



[Hackman et al.|2012| ). Early DI images corresponding to the pe- 
riod between 1994 and 2002 suggested persistent presence of a 
pair of active longitudes and showed major changes in the sur- 
face structure on a time-scale of less than a year (Berdyugina[ 
et al.|1998a|p^99l[Lindborg et al.^201 1 ). More recent DI maps 



Roettenbacher et al. 201 1 , and references therein), aiming to ex- 



plore long-term activity cycles and, in particular, to investigate 



covering the years from 2004 to 2010 revealed the star entering 
a low-activity state characterised by a more random distribution 
of cool spots ( Hackn mi et al. 1201 2]). 

The underlying cause of the remarkable surface activity of 
II Peg - the d ynamo-generated ma gnetic field - was first detected 
in this star by [Donati et al.| ( p^92l ) with four circular polarisation 
observations of a few magnetically sensitive lines. The presence 
of the Stokes V signatures in spectral lines was subsequently 
confirmed by [Donati et al.[ ( [1997| ) using a multi-line polarimetric 
analysis, but no systematic phase-resolved investigation of the 
magnetic field topology of II Peg has ever been undertaken. 

Since 2004 we have been monitoring the magnetic field in 

11 Peg using the SOFIN spectropolarimeter at the Nordic Opti- 
cal Telescope. We have acquired a unique collection of high- 
resolution Stokes / and V spectra covering ^^5.5 years or al- 
most 300 stellar rotations. In this paper we present a compre- 
hensive analysis of these polarisation data, focusing on the self- 
consistent ZDI mapping of stellar magnetic field topology. The 
accompanying study by Hackman et al.[ ( [2QT2l ) obtained temper- 
ature maps from the same SOFIN Stokes / spectra. Preliminary 
attempts to map the magnetic field geometry of II Peg, using 
about 10% of the spectropolarimetric data analysed here, were 
presented fi rst by [Carroll et al. (2007) , and then by Carroll et al. 
( [2009b| ) and [Kochukhov et al.[ ( [200'9f . 

This paper is organised as follows. In Sect. [2] we describe 
the acquisition and reduction of the spectropolarimetric observa- 
tions of II Peg. Detection of the magnetic signatures in spectral 
lines with the help of a multi-line analysis is presented in Sect. [3] 
The methodology of the self-consistent ZDI and the choice of 
stellar parameters required for mapping is described in Sect. [4] 
Results of the magnetic and temperature inversions of II Peg are 
presented and analysed in Sect. [5] The outcome of our investi- 
gation is discussed in the context of previous observational and 
theoretical studies in Sect. O 

2, Spectropolarimetric observations 

The spectropolarimetric observations of II Peg analysed here 
were carried out during the period from Jul. 2004 to Jan. 2010 
using the SOFIN echelle spectrograph ( [Tuominen et al.[[T999[ ) 
at the 2.56-m Nordic Optical Telescope. The spectrograph, 
which is mounted at the Cassegrain focus and is equipped with 
a 2048 X 2048 pixel CCD detector, was configured to use its 
second camera yielding resolving power of 70000. For ZDI 
analysis we used 12 echelle orders, each covering 40-50 A in the 
wavelength region between 4600 and 6135 A. The data for II Peg 
were obtained during 12 individual epochs, for which from 3 to 

12 observations were recorded over the time span ranging from 
3 to 14 nights. The spectra have typical S/N ratio of 200-300. 

The circular polarisation observations were obtained with the 
Zeeman analyser, consisting of a calcite plate used as a beam 
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splitter and an achromatic rotating quarter-wave plate. At least 
two exposures with the quarter-wave retarder angles separated 
by 90° is required to obtain the Stokes V spectrum. Rotation of 
the quarter-wave plate has the effect of exchanging positions of 
the right- and left-hand circularly polarised beams on the detec- 
tor. The beam exchange procedure ( |Semel et al.|1993[ ) facilitates 
an accurate polarisation analysis because possible instrumental 
artefacts change sign when the quarter- wave plate is rotated and 
then cancel out when all sub-exposures are combined. To ac- 
cumulate sufficient signal-to-noise ratio, a sequence of 2 or 3 
double exposures was obtained in this way. The length of indi- 
vidual sub-exposures varied between 15 and 25 min, depending 
on seeing and weather conditions. 

The data were reduced with the help of the 4A software 
package plyin||200Q| ). Specific details of the SOFIN polarime- 
ter design and corresponding data reduction methods are given 
by |Ilyin| ( |2012| ). The spectral processing included standard re- 
duction steps, such as bias subtraction, flat field correction, re- 
moval of the scattered light, and optimal extraction of the spec- 
tra. Wavelength calibration used ThAr exposures obtained be- 
fore and after each single exposure to account for environmen- 
tal variations and gravitational bending of a Cassegrain mounted 
spectrograph by means of a global fit of the two ThAr wave- 
length solutions versus time. 

Spectropolarimetric observations of strongly magnetic Ap 
stars are frequently performed with SOFIN in a configuration 
similar to the one used for II Peg. These measurements agree 
closely with the results obtained for the same stars at other tele- 
scopes (e.g. |Ryabchikova et al.|2007]|Ilyin|20T2] ). Repeated ob- 
servations of a very slowly rotating Ap star y Equ reveals no 
systematic differences in the Stokes V profiles during the entire 
period of our observations of II Peg. This confirms the robust- 
ness and accuracy of the employed instrument calibration and 
data reduction methods. 

The log of all 88 SOFIN Stokes V observations of II Peg 
is given in Table [T] The rotation of the primary component of 
II Peg is synchronised with the orbital motion. Therefore, we 
calculated rotational phases using the orbital ephemeris from 
[Berdyug inaetaL] ( [1998b| 

T = 2449582.9268 -h 6.724333 x E (1) 

which refers to the time of orbital conjunction. The orbital solu- 
tion from the same study was applied to correct the radial veloc- 
ity variation caused by the orbital motion. 

Our 12 sets of spectropolarimetric observations of II Peg 
comprise 3 to 12 distinct rotational phases, whereas -10 evenly 
distributed phases are needed for an optimal ZDI reconstruction 
( [Kochukhov & Piskunov"2002). Assuming that each observa- 
tion provides a coverage of 10% of the rotational period, we 
quantified the phase coverage of each data set as a fraction / of 
the full rotational cycle (see Table [T]). Our SOFIN observations 
yield from / = 30% (epoch 2005.0) to / = 87% (epoch 2009.7). 
The three epochs with the best phase-coverage are 2004.6 (75%), 
2007.6 (81%), and 2009.7 (87%). 

The quality of ZDI mapping is reduced for the data sets with 
an insufficient phase-coverage. Yet, as demonstrated by numeri- 
cal experiments ( Donati & Brown 1997 ; Kochukhov & Pis kunov] 
[2002 ) and ZDI studies bas ed on observatio ns with poor phase- 
coverage (Donati, 1 999^ ^Hussain et al.|2009| ), it is still possible to 
extract useful information about the stellar surface features from 
just a few spectra and from observations covering only half of 
the stellar rotation cycle. Aiming to maintain consistency in our 
long-term ZDI analysis of II Peg, we performed magnetic inver- 
sions for all 12 epochs. At the same time, we kept in mind the 



diflference in data quality during the assessment of the inversion 
results. 



3. Multi-line polarisation analysis 

3.1. Least-squares deconvolution 

Even the most active late-type stars exhibit a relatively low am- 
plitude of the circular polarisation signal in spectral lines, rarely 
exceeding 1% of the Stokes V continuum intensity. Given the 
moderate quality of the SOFIN spectropolarimetric observations 
of II Peg, we were able to detect Stokes V signatures at the 2- 
3(T confidence level in only a few of the strongest spectral lines 
(Fig.[T]). This makes it challenging to model the magnetic field 
of II Peg using a direct analysis of individual line profiles. A 
widely-used approach to overcome this difficulty is to employ a 
multi-line technique, combining information from hundreds or 
thousands of individual metal lines. 

In this study of II Peg we applied the least- squares decon- 
volution (LSD) technique f Donati et al.| 1997 ) us ing the code 
and methodology described by [Kochukhov et aL] ( |2010| ). The 
LSD technique extracts information from all available lines by 
assuming that the Stokes / and V spectra can be represented by 
a superposition of corresponding scaled mean profiles. The scal- 
ing factors, established under the weak-line and weak-field ap- 
proximations, are equal to the central line depth d for Stokes / 
and to the product of the line depth, laboratory wavelength of 
the transition A, and its eflTective Lande factor z for Stokes V, 
respectively. A linear superposition of scaled profiles is mathe- 
matically equivalent to a convolution of the average profile with 
a line mask. This simple stellar spectrum model can be inverted, 
obtaining a high-q uality mean profile for a given line mask and 
observational data. [Kochukhov et al.| ( |2010f showed that, for the 
magnetic fields below ^ 2 kG, the LSD Stokes V profiles derived 
in this way can be interpreted as a real Zeeman triplet line with 
an average Lande factor. 

The atomic line data required for the application of LSD 
to observations of II Peg were obtained from the vald database 
( |Kupka et aL]|1999| ). We extracted a total of 1580 spectral 
lines with the central depth larger than 10% of the contin- 
uum. The line intensities were calculated for a marcs model 
atmosphere ( [Gustafsson et al.|[2008| ) with the eflTective temper- 
ature Teff = 4750 K, surface gravity log g = 3.5, and metallicity 
[M/H] = -0.25. The mean parameters of the resulting line mask 
are Aq = 5067 A, zo = 1.21, and do = 0.46. The same parame- 
ters were used for normalising the Stokes / and V LSD weights. 
The LSD profiles were calculated with a 1.2 kms"^ step for the 
velocity range ± 60 km s"^ 

The LSD analysis of the SOFIN observations of II Peg al- 
lowed us to achieve slS/N gain of 30-40, obtaining a definite de- 
tection (false alarm probability < 10"^) of the circular polarisa- 
tion signatures for all but one Stokes V spectrum (phase 0.035 for 
epoch 2008.7). The resulting mean Stokes V profiles are compat- 
ible with the marginal polarisation seen in individual lines (see 
Fig.[T]), but have a much higher quality, making them suitable for 
detailed modelling. 



3.2. Mean longitudinal magnetic field 

In Table [T] we report the maximum amplitude of the LSD Stokes 

V profiles and the mean longitudinal magnetic field (^z) deter- 
mined from their first moment. We found typical LSD Stokes 

V amplitudes of the order of 0.1% and (^z) in the range 30- 
250 G, determined with the precision of 5 G. These results are 
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Fig. 1. Comparison of the Stokes V profiles of the Cai 6122.2 A 
line (left panel) with the scaled LSD V profiles (right panel) at epoch 
2004.6. The spectra for diff'erent rotational phases are shifted vertically. 
The 1-cr error bars are indicated on the left side of each panel. 



compatible with the previous longitudinal field measurements by 



^gt]([T98a), which had an accuracy of 100-160 G and did not 
yield a convincing detection of magnetic field in II Peg. 

Despite the ubiquitous definite detections of the circular po- 
larisation signatures, 16 out of 88 LSD profiles of II Peg do not 
exhibit a significant mean longitudinal magnetic field. A com- 
parison of (5z) measurements with the maximum LSD Stokes 
V amplitude (Fig. [2]) shows no obvious correlation. This in- 
dicates that, with a few exceptions (e.g. phase 0.531 at epoch 
2004.6 illustrated in Fig.[T]), the Stokes V profiles of II Peg are 
generally complex and cannot be fully characterised by their 
first-order moment. Nevertheless, analysis of the longitudinal 
magnetic field measurements does reveal an intriguing long-term 
behaviour. As illustrated in Fig. [3] the range of the longitudi- 
nal field variation with rotational phase systematically decreased 
from 2004 to 2010. Furthermore, it appears that the phase- 
averaged value of the longitudinal magnetic field was signifi- 
cantly deviating from during the period from the beginning of 
our monitoring until the epoch 2006.9. This suggests the pres- 
ence of large areas of positive magnetic polarity on the visible 
hemisphere of the star before that epoch and a more symmetric 
distribution of the field orientation at later epochs. 

4. Self-consistent magnetic and temperature 
mapping 

4.1. Zeeman Doppler imaging code InversIS 

Reconstruction of the distribution of magnetic field and tem- 
perature on the surface of II Peg was carried out with the help 
of our new ZDI code Invers13. This inversion so ftware repre- 
sents a further development of the InversIO code ( [Piskunov~& 



|Kochukho v'2002>'Kochuk hov & Piskunov|2002| ), previously ap- 
plied for mapping magnetic geometries and chemical spots on 
Ap stars ([Kochukhov et"aLl[2004l [Kochukhov "&Wa de 2010; 
|Luftinger et al.|2010| ). Both codes incorporate full treatment of 
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Fig. 3. Range of longitudinal field variation with rotational phase 
for diff'erent epochs of II Peg observations. The line thickness is pro- 
portional to the phase coverage /. The squares show longitudinal field 
averaged over rotational phase. 



the polarised radiative transfer using a modified Diagonal El- 
ement Lambda-Operator (DEL Q) algorithm (Rees et aL][T989) 
Piskunov & Kochukhov||2002| ). All four Stokes parameters are 
calculated simultaneously and fully self-consistently, taking into 
account Zeeman broadening and splitting of spectral lines in the 
intensity spectra on the one hand, and attenuation of the po- 
larised profiles due to temperature spots on the other hand. Some 
aspects of this approach have already been incorporate d in the 
ZDI of cool stars in the preliminary studies of II Peg ( jCarroll| 
etal.|200 7i'. 200%; K ochukhov et al. 2009 ). However, the major- 
ityof recent applic ations of ZDI (e.g. Skelly et al. 2010; Waite| 
et al.|2011y ^Marsden et al.|20lT] ) still systematically neglect ef- 
fects of temperature spots in the magnetic inversions. As demon- 
strated by the numerical tests of Rosen & Kochukhov (20 12), the 
lack of self-consistency in temperature and magnetic mapping 
can lead to severe artefacts if cool spots coincide with major 
concentrations of the magnetic flux. 
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Calculation of the local Stokes parameter profiles by In- 
VERsl3 is based on a prescribed line list of atomic and molecu- 
lar lines, normally obtained from the vald and marcs databases, 
and a grid of model atmospheres. In the present study of II Peg 
we employed a grid consisting of 18 marcs model atmospheres 
covering a T^^ range of 3000-5750 K with a 100-250 K step 
in temperature. The Stokes IQUV profiles and continuum in- 
tensities Ic are calculated for the temperature of a given stellar 
surface element by quadratic interpolation between three sets of 
model spectra corresponding to the nearest points in the model 
atmosphere grid. This allows an accurate semi- analytical com- 
putation of the derivatives with respect to temperature. On the 
other hand, the derivatives with respect to the three magnetic 
field vector components are evaluated numerically, using a sim- 
ple one-sided diflTerence scheme. 

The local profiles are convolved with a Gaussian function 
to take into account instrumental and radial-tangential macrotur- 
bulent broadening, Doppler shifted, and summed for each rota- 
tional phase on the wavelength grid of observed spectra. The 
resulting Stokes parameter profiles are normalised by the phase- 
dependent continuum flux and compared with observations. In 
calculating the goodness of the fit we approximately balance the 
contributions of each Stokes parameter by scaling the respective 
chi- square terms by the inverse of the mean amplitude of the 
corresponding Stokes profiles. Test inversions by Kochukhov & 
[Piskunov ( 2002 ) and Rosen & Kochukhov ( 2012 ) verified simul- 
taneous reconstruction of the magnetic and starspot maps using 
this weighting scheme. This approach is also routinely used in 
the reconstruction of Ap-star magnetic fields from four Stokes 
parameter observations ("Kochukhov et al. 2004[ [Kochukhov & 



[Wade 2010 ) . In the study of II Peg the contributions of the Stokes 
/ and V spectra were weighted as 1:30. Inversion results are not 
sensitive to the exact choice of the relative weighting of Stokes 
/ and V. 

The iterative adjustment of the surface maps for match- 
ing available observations is accomplished by means of the 
Levenberg-Marquardt optimisation algorithm, which enables 
convergence in typically 10-20 iterations starting from a ho- 
mogeneous temperature distribution and zero magnetic field. 
The Invers13 code is optimised for execution on massively- 
parallel computers using MPI libraries. We refer the reader to 
[Kochukhov et al.] ( [201 2| ) for further details on the numerical 
methods and computational techniques employed in Invers13. 

As shown numerically by [Donati & Brown ( 1997 ) and 
[Piskunov & Kochukhov] p002| ), and analytically by [Piskunov 



ber m (see Eqs. 2-8 in [ Donati et al.[[2006| ). In particular, the 
radial field is determined entirely by the poloidal contribution 
(coeflftcient a), whereas both poloidal and toroidal components 
(coeflftcients JS and y, respectively) are contributing to the merid- 
ional and azimuthal field. In this description of the stellar mag- 
netic field topology through a spherical harmonic expansion, the 
role of regularisation is played by a penalty function equal to 
the sum of squares (magnetic energies) of the harmonic coeflft- 
cients weighted by This prevents the code from introducing 
high-order modes not justified by the observational data. 

Expansion in spherical harmonics has the advantage of au- 
tomatically satisfying the divergence-free condition for a mag- 
netic field. It also allows us to conveniently characterise the rela- 
tive contributions of diflTerent components to the stellar magnetic 
field topology and to study the evolution of these contributions 
with time. Therefore, for the magnetic mapping of II Peg we 
adopted the spherical harmonic approach, choosing ^^ax = 10 
which corresponds to a total of 360 independent magnetic field 
parameters. 

Similar to other quantities obtained by solving a regularised 
ill-posed problem, the spherical harmonic coefl&cients cannot be 
attributed formal error bars because, in addition to the observa- 
tional data, the individual modes are constrained by the penalty 
function. Therefore, their formal error bars are strongly depen- 
dent on the choice of regularisation parameter, especially for 
the high-^ modes. This is why here, as well as in all previous 
ZDI studies employing spherical harmonic expansion (e.g. [Do-[ 
nati et al.|20Q6l|Skelly et al.|2010||F^es et al.|2012| ), no formal 



uncertainty estimates are derived for the harmonic coeflftcients. 

As inferred from the previous numerical tests of ZDI inver- 
sions ([Donati & Brown [T997[ [Kochuk hov & Piskunov] [20021 



Rosen & Kochu khov 2012]), reconstruction of the magnetic field 
distributions from the Stokes / and V data alone suffers from 
some ambiguities. The most important problem is a cross-talk 
between the radial and meridional field components at low lati- 
tudes (i.e. for the regions below the stellar equator for the rota- 
tional geometry of II Peg). On the other hand, the harmonic reg- 
ularisation or parameterisation of the magnetic field maps partly 
alleviates this problem (Donati 2001; Ko chukhov & Piskunov] 
[2002) ), because in this case the inversion code is not allowed to 
vary the three magnetic vector components arbitrarily. Instead, 
the meridional and azimuthal components are coupled via the JS 
and y harmonic coefl&cients which avoids the cross-talk between 
the radial and meridional fields, at least for the low-^ values. 



(|2005| ZDI with oHlTstokes / and V spectra is an intrinsi- ^ 2. Reconstructing temperature spots from LSD profiles 



cally ill-posed problem, requiring the use of a regularisation or 
penalty function to reach a stable and unique solution. The re- 
construction of temperature in Invers13 is regularised using the 
Tikhonov method, very similar to many previous DI studies of 
cool active stars (e.g. (Piskunov & Rice|1993| ) and our recent tem- 
perature mapping of II Peg with the same SOFIN data ( [Hackman 



[et al. 2012 ). In this approach the code seeks a surface distribution 
with a minimum contrast between adjacent surface elements. 

For the magnetic field reconstruction, Invers13 offers a 
choice of either directly mapping the three magnetic vector com- 
ponents and applying the Tikhonov regularisation individually to 
the radial, meridional, and azimu thal magnetic component maps 
( Piskunov & Kochukhov|[2002[), or expa nding the field into a 
spherical harmonic series ( [Donati et al.|2 006 ). In the latter for- 
mulation of the magnetic inversion problem, the free parameters 
optimised by the code are the spherical harmonic coefl&cients 
0!{m^ Ptm^ and 7^^, giving the amplitudes of the poloidal and 
toroidal terms for a given angular degree i and azimuthal num- 



Most previous ZDI modelling of the Stokes / and V LSD profiles 
of late-type active stars relied on simple analytical formulas, of- 
ten a Gaussian approximation, to represent the local intensity 
and polarisation profiles. These coarse line-profile models have 
either entirely neglected dependence of the local line intensity 
on temperature (e.g. [Petit et al.| |2004 ) or used disk-integrated 
spectra of suitably chosen cool and hot slowly rotating template 
stars to represent local contribution of the spot and photosphere, 
respectively (e.g.lD onati et al.|20Q3j ). 

Here we develop a more sophisticated approach to inter- 
pret the LSD profile time series. We aim to describe the LSD 
Stokes / and V spectra with self-consistent polarised radiative 
transfer calculations, based on realistic model atmospheres, for 
a fictitious spectral line with average atomic parameters. Com- 
prehensive investigation of the properties of LSD profiles by 
[Kochukhov et al.[p010| showed that this approach yields an ac- 
curate description of the mean circular polarisation line shapes 
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Fig. 4. Comparison of temperature mapping from a real line and 
from LSD profiles, a) True temperature distribution, b) DI map re- 
constructed from the simulated observations of a real spectral line, c) 
Temperature reconstruction from the simulated LSD profiles assuming 
mean-line parameters. 



for the magnetic field strengths below ^ 2 kG. For stronger fields, 
details of the LSD Stokes V profile shape cannot be reproduced, 
although the LSD profile moments still yield correct longitu- 
dinal magnetic field. According to previous studies, magnetic 
fields found on the surfaces of RS CVn- and BY Dra-type active 
stars are safely within this limit. On the other hand, analyses 



of Ap/Bp stars hosting much stronger fields (e.g. Wade et al 
2000} [Silvester et aL][MT2| ) proves that the LSD line-averaging 



method itself does not systematically underestimate magnetic 
field strength even beyond the weak-field regime. 

At the same time, [Kochukhov et al.| PQ1Q| ) found that a re- 
sponse of the LSD Stokes / profile to the variation of chemical 
composition differs from that of any real spectral line, which 
prevents an accurate mapping of the chemical spots in Ap stars 
using LSD profiles ( [Folsom et al.||2QQ8|| ). One might suspect 
that reconstruction of the temperature inhomogeneities from the 
LSD spectra is similarly limited. It is clear that the tempera- 
ture sensitivity of an average profile, composed of an essentially 
random mixture of thousands of lines with diverse temperature 
responses, cannot be known as well as the one for a set of judi- 
ciously chosen individual spectral lines with accurate atomic pa- 
rameters. On the other hand, typical Stokes / atomic line profile 
variations in active stars are dominated by the continuum bright- 
ness effect, which influences all lines in a sim ilar way ( |Vogt &| 
Penrod|1983]|Unruh & Collier Cameron|p^ 5). Adding to this 
a greatly enhanced S/N of the LSD spectra, one might hope to 
achieve a reasonable reconstruction of temperature spots treat- 
ing the LSD intensity profiles as a real spectral line. To clarify 
whether this is indeed possible, we have carried out numerical 
tests of temperature inversions based on simulated LSD spectra. 

First, we established mean-line parameters appropriate for 
the LSD line mask applied to the SOFIN observations of II Peg. 
Among the 1580 spectral lines included in this mask, Fei is the 



most common ion (about 40% of all lines). Thus, we adopted the 
Fe I identification for a fictitious line representing LSD profiles 
and obtained the preliminary set of atomic parameters necessary 
for spectrum synthesis by combining the average excitation en- 
ergy of the lower level, oscillator strength, and broadening con- 
stants for this ion with the mean wavelength and Lande factor 
(5067 A and 1.21, respectively) for the entire line list. Then 
we adjusted the oscillator strength and van der Waals damping 
constant by fitting the synthetic LSD profiles derived from a the- 
oretical spectrum covering the same wavelength interval as real 
observations. These calculations adopted the stellar atmospheric 
parameters discussed below (Sect. |4.3| ) and were performed for 
the Ve sin / of II Peg and for the non-rotating star with the same 
parameters. In both cases we found that single-line calculations 
provide an excellent fit to the LSD profiles obtained from the 
synthetic spectra. 

At the next step, we simulated spectra of a star with tem- 
perature inhomogeneities, using a photospheric temperature 
of 4750 K and a spot temperature of 3700 K. The Stokes 
/ spectra for the entire SOFIN wavelength range were com- 
puted for 10 equidistant rotational phases, assuming / = 60°, 
Ve sin/ = 23 kms"\ and no magnetic field. A surface tempera- 
ture map consisting of four circular spots located at diflTerent lat- 
itudes (see Fig. |4^) was employed in these calculations. These 
theoretical spectra were processed with our LSD code in exactly 
the same way as real stellar observations. The resulting set of 
synthetic LSD profiles was interpreted with Invers13 using the 
line parameters established above. For comparison, we also re- 
constructed a temperature map from the direct forward calcula- 
tions for the same spectral line. 

A comparison of the DI maps inferred from the simulated 
data with the true temperature distribution is presented in Fig. [4] 
The reconstructions from the simulated observations of the real 
spectral line and from the simulated LSD profiles are generally 
successful in recovering the properties of cool spots. However, 
the temperature DI using LSD profiles shows a tendency to pro- 
duce spurious hot surface features in which the temperature is 
overestimated by 200-300 K. We emphasise that this problem 
appears entirely due to a loss of information about real temper- 
ature sensitivity of the local profiles of individual spectral lines. 
Fortunately, this bias does not significantly interfere with the re- 
construction of large low-temperature spots, although their con- 
trast seems to be systematically underestimated. Thus, we con- 
clude that the temperature DI based on a single-line radiative 
transfer approximation of the LSD profiles is a viable approach 
to mapping the distribution of cool spots on the surfaces of active 
stars. At the same time, the reliability of hot features recovered 
with this method is questionable. 

Limitations of LSD-based temperature inversions partly ex- 
plain why some of the temperature maps presented her e diflTer 
from those in our previous analysis ( [Hackman et al. P2012 ). The 
diflferences are mainly caused by the use here of LSD profiles, 
while our previous analysis was based on individual lines. While 
individual lines will give a more accurate result, this approach 
suflTers from a lower S/N. At the same time, neglect of the Zee- 
man broadening and intensification in earlier temperature inver- 
sions may also contribute to the discrepancy with the present 
results. In general, both of these eflfects will influence the spot 
temperatures and latitudes. 
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4.3. Stellar parameters 

DiflPerent sets of atmospheric parameters have been deduced in 
previous studies of II Peg. The most thorough model atmo- 
sphere analyses were those by IBerdyugina et al.| ( |1998 b|) and 
Qttmann et aLl ( |1998| ). The first paper suggested T^^ = 4600 K, 
log ^ = 3.2, and [M/H] = -0.4, whereas the second paper has in- 
ferred Teff = 4800 K, log^ = 3.65, and [Fe/H] = -0.24. A com- 
parison of the average SOFIN observations with the theoretical 
spectra computed using the synth3 code ( Kochukhov 2007 ) and 
MARCS model atmospheres supports the latter set of stellar pa- 
rameters. Therefore, we adopted a marcs model atmosphere 
with Teff = 4750 K, log g = 3.5, and [M/H] = -0.25 to represent 
an unspotted star. This temperature was also used as a start- 
ing guess for DI. All spectrum synthesis calculations assumed a 
microturbulent velocity of = 2 kms"^ and a radial-tangential 
macroturbulent broadening of ft = 4 kms"^ These parameters 
are compatible with the values determined in previous studies. 
We note that the synthetic Fe i line approximating the LSD pro- 
files of II Peg is weak (W;i = 40 mA) and, therefore, relatively 
insensitive to the choice of 

Based o n previous studies (Berdyugi na et al.|1998b[ |Frasca 
et al.||2008 ]), we adopted an inclination of the stellar rotational 
axis / = 60°. Using the three data sets with the best phase- 
coverage (epochs 2004.6, 2007.6, and 2009.7), we found that the 
best fit to the Stokes / LSD profiles is achieved for Ve sin / = 23 ± 
0.5 km s"^ This value of the projected rotational velocity is con- 
sistent with Ve sin / = 22.6+0.5 km s"^ determined by Berdyugina^ 



et al.|(^99 8b) and used in the recent DI studies by |Lindborg et al. 
( 12011 



and Hackmanet al.|(2012). 



A latitude-dependent diflferential rotation can be incorpo- 
rated in the ZDI with Invers13. This effect was ignored by pre- 
vious spectroscopic DI investigations of II Peg and was theoreti- 
cally predicted to be negligi ble in the rapidly rotating stars with 
deep convective envelopes ( [Kitchatinov & Rudiger| 1999] ). On 
the other hand, photometric analyses occasi onally claime d de- 
tection of the diflTerential rotation in II Peg ( Henry et al.|p^95| 
[Siwak et al.,^2010( |Roettenbacher et al.||2011| ), yet giving con- 
tradictory results, mainly because of the diflSculty of constrain- 
ing the spot latitudes in light curve modelling. However, even 
the largest value of the diflTerential rotation coefl&cient reported 
in the literature k = 0.0245 ( |Siwak et al.||2010| ) is substantially 
smaller than the solar value of ^ = 0.19 and corresponds to an 
insignificant surface shearing for the time span of most of our 
SOFIN data sets. Therefore, ZDI analysis of II Peg was carried 
out assuming no diflTerential rotation. 



5. Magnetic field topology of II Peg 

Results of the simultaneous reconstruction of magnetic field 
topology and temperature distribution for II Peg are illustrated 
for the four data sets with best phase-coverage in Figs. [5] and 
[6] and for the remaining epochs in the online Figs. |7}|9| For 
each data set we present four rectangular maps, showing the dis- 
tribution of temperature as well as the radial, meridional, and 
azimuthal field components. A comparison between observed 
and computed LSD Stokes / and V profiles is illustrated next to 
the corresponding surface images. The model profiles generally 
achieve an adequate fit to observations, with median mean devi- 
ations of 0.15% and 0.01% for the Stokes / and V LSD profiles, 
respectively. 

The maximum local magnetic field intensity inferred by our 
ZDI code reaches 1 kG for 2004.6 data set, but is typically on 
the order of 400-500 G for the later magnetic images. Judging 



from these results, the field strength on the surface of II Peg was 
steadily decreasing from our first observations to later epochs. 
The field topology appears to be rather complex in 2004-2005, 
with many localised magnetic spots in all three magnetic vector 
component maps. As the field became weaker, this configura- 
tion evolved into a simpler topology, often dominated by large 
regions of common field polarity (e.g. epoch 2007.6). 

The temperature maps of II Peg reconstructed from LSD pro- 
files show cool spots with an 800-1000 K contrast relative to 
the photospheric T = 4750 K. Occasionally, the inversion code 
reconstructed spots which are a few hundred K hotter than the 
photosphere. However, the numerical experiments discussed in 
Sect. |4.2| suggest that these features are not reliable. 

Inversions were carried out for all epochs, including the 
three data sets (2005.0, 2005.9, and 2008.7) with a poor phase- 
coverage. Temperature maps for these three data sets are dom- 
inated by a strong axisymmetric component, which should be 
regarded as spurious. 

5.1. Average field parameters 

The two independent ZDI reconstructions obtained for the 
epochs 2006.7 and 2006.9 allow us to investigate short-term 
changes of the surface distributions. The temperature maps in- 
ferred from the observations separated by three months are very 
similar. Both images are dominated by a large low-latitude cool 
spot at 180° longitude and an extended high-latitude spot group 
in the longitude interval 300° to 60°. The azimuthal field maps, 
representing the dominant magnetic component for both epochs, 
are also fairly similar. At the same time, distributions of the 
weaker radial and meridional field components exhibit signifi- 
cant diflTerences, suggesting field evolution. 

Next, to detect possible systematic changes either in the 
phases of the surface magnetic field distribution or in the total 
magnetic field strength, we calculated diflTerent types of aver- 
ages of the ZDI maps. First, we averaged over the visible lati- 
tude range in each map, to obtain magnetic field profiles that de- 
pend on longitude (or phase) only. The profiles shown in Fig.fTT 
are averaged over the latitude interval from -60° to 90° for the 
diflTerent magnetic field components as function of time. This 
procedure can reveal azimuthal dynamo waves, that were visible 
in th e surface temperature maps of this object during 1994-2002 
(see |Lindborg et alT||2011j ). There is no evidence for this type 
of dynamo wave from the ZDI maps. This agrees with the re- 
sults of |Hackman et al.| ( |2012| ), who reported the disappearance 
of the clear drift pattern in the surface temperature maps during 
these years. The radial field plot, on the other hand, reveals a 
rather abrupt appearance and disappearance of spots of opposite 
polarities at a certain phase which are, however, irregular with 
time. 

Finally, we calculated the root-mean-square values of the 
magnetic field over each ZDI map, characterising the overall 
magnetic field strength (shown in Fig. [12]). Interestingly, all the 
magnetic field components are of comparable strength, immedi- 
ately hinting, from basic dynamo theory, towards an dynamo 
operational in the object - the presence of diflTerential rotation 
would lead to the eflftcient shearing of the poloidal field into 
a toroidal field, in which case the azimuthal component would 
be observed to dominate over the other components. The ra- 
dial and meridional components are clearly decreasing monoton- 
ically with time, reaching a minimum at around the year 2009. 
After that, the radial field possibly starts rising again. The az- 
imuthal field shows a somewhat diflferent trend in time: at first it 
is somewhat weaker than the other components, slightly increas- 
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Fig. 5. Magnetic field topology and temperature distribution of II Peg at epochs 2004.6 and 2006.7. For each epoch the two columns on the 
left compare the observed Stokes / and V LSD profiles (symbols) and the theoretical spectrum synthesis fit (solid lines). Spectra corresponding 
to diff'erent rotational phases are shifted vertically. The scale is given in % for Stokes V and in units of the continuum intensity for Stokes /. 
The four rectangular maps represent surface distributions of the radial, meridional, and azimuthal magnetic field components and corresponding 
temperature reconstructed for each epoch. The contours are plotted with a step of 0.2 kG over the magnetic maps and with a step of 200 K over 
the temperature distributions. The thick contour lines indicate where the field changes sign in the magnetic maps. The vertical bars above each 
rectangular panel indicate rotational phases of individual observations. The rotational phase runs from right to left on all maps. 



ing during 2004-2007, after which it also shows a decreasing 
trend. Similar to the radial component, it seems to start rising 
again after 2010. Therefore, judging from the magnetic field 
strength, it seems plausible that the magnetic activity level of 
the star has been declining during the epoch 2004-2009, while 
the signs of rising activity can be seen at least in the rms field 
strength for later epochs. 

5.2. Evolution of harmonic field components 

Magnetic field maps reconstructed by our ZDI code are param- 
eterised in terms of the spherical harmonic coeflftcients corre- 



sponding to the poloidal and toroidal field components. In ad- 
dition to the analysis of 2D maps, this representation provides 
another convenient possibility to characterise the field topology 
and its long-term evolution. 

The sum of the spherical harmonic coefiftcients squared is 
proportional to the total energy contained in the stellar magnetic 
field. Figure [TSh illustrates how this parameter has changed for 
II Peg between 2004 and 2010. Consistently with the results dis- 
cussed above, we find a noticeable decrease of the total field 
energy in the period between 2004 and 2008, with a possible re- 
versal of this trend afterwards. One can note that the total mag- 
netic energy recovered from the data sets containing only a few 
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Fig. 6. Same as Fig.|5]for epochs 2007.6 and 2009.7. 



spectra (e.g. 2005.0, 2005.9) is systematically underestimated [Fares et aL]|2009| ). With this definition, non-axisymmetric field 



compared to the trend defined by other maps. It is reassuring 
that the analysis of these data sets does not result in spurious 
strong magnetic field features. 

The time dependence of the relative contributions of the 
poloidal and toroidal field components is illustrated in Fig. \T3\ ). 
It appears that our observations of II Peg reveal a cyclic change 
of the field topology on the time scale of a few years. Ignoring 
magnetic maps corresponding to epochs with a partial phase- 
coverage, one can conclude that II Peg exhibited a predominantly 
poloidal field before 2007.6 and a mainly toroidal field after- 
wards. However, the difference between the energies of the two 
components is never much larger than 20-30%. 

Finally, Fig. [13]: assesses relative contribution of the axisym- 
metric and non-axisymmetric field components as a function of 
time. Here we define axisymmetric harmonic components as 
those with m < ijl and non-axisymmetric ones as m > £12 (e.g. 



dominates all the time, except for the last two observing epochs. 



As mentioned in Sect. 4.1 formal error bars cannot provide 
a realistic estimate of uncertainties of the spherical harmonic co- 
eflftcients recovered in a regularised least-squares problem. In- 
stead, the scatter of points corresponding to close observational 
epochs gives an idea of the uncertainties. We can see that in 
Figs. [12 and[T3]the points inferred from the poor phase-coverage 
data sets often deviate significantly from the general trends. On 
the other hand, there are only a few cases when results obtained 
from the good phase-coverage data exhibit abrupt changes. This 
suggests that the trends examined in this and the previous sec- 
tions are real and are not dominated by random inversion errors. 
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Fig. 11. Magnetic field components (top: radial field, middle: merid- 
ional field, bottom: azimuthal field) averaged over the latitude interval 
-60° to 90° from each ZDI map. Each stripe represents the average 
field from one ZDI map, extended over the time axis to help the visu- 
alisation. The stripes are organised as function of time, according to 
their observational epoch, running from top to bottom. Dark colours 
represent negative polarities, bright colours positive ones; in all figures 
a linear colour table with minimum of -300 G and maximum of 300 G 
has been used. 



5.3. Extended magnetospheric structure 

The results of ZDI calculations are commonly used to investi- 



gate an extended stellar magnetospheric structure (Donati et al. 
20Q8l|Jardine"etal . 2008; G regory et al.]2008l [Fares et al.|2012| ). 



The knowledge of the field topology above the stellar surface and 
in the immediate circumstellar environment allows photospheric 
magnetic field measurements to be connected with the studies of 
stellar coronas, X-ray emission, and prominences and the inter- 



action between the mass loss and magnetic field to be investi- 
gated. To determine the structure of the circumstellar magnetic 
field, one can use the potential field source surface (PFSS) ex- 
trapolation method devel oped for the solar magnetic field by |van| 
[Ballegooijen et al.|( 1998 ) and adapted for stellar magnetic fields 
by |Jardine et al.| ( |2002| ). In this method the extended stellar mag- 
netic field is represented as the gradient of a scalar potential ex- 
panded in a spherical harmonic series. The boundary condition 
for the radial field component at the stellar surface is provided 
by the empirical magnetic field maps. The outer boundary con- 
dition is given by a source surface of radius Rs beyond which the 
field is assumed to be purely radial. 

We reconstructed the extended magnetospheric structure of 
II Peg with the help of an independently developed PFSS code. 
This software was applied to all our ZDI maps. The source sur- 
face is placed at = 37?^, which is plausible given the mean 
value of the solar source surface radius of 2.5Rq. Previous po- 
tential field extrapolation st udies adopted similar R^ values for 
other cool active stars (e.g. |Jardine et al.||2002] |Hussain et al 
[2002] ). 

The magnetospheric structure of II Peg predicted by our ZDI 
maps is illustrated in Fig.[T4]for the four epochs with best phase- 
coverage. The stellar magnetosphere is shown at four distinct 
rotational phases, with diflTerent colours highlighting open and 
closed magnetic field lines. The evolution of the large-scale 
field is evident from this figure. Since the contribution of the 
small-scale (and hence less reliably reconstructed) field struc- 
tures decays more rapidly with radius, this potential field extrap- 
olation essentially presents a distilled view of the radial compo- 
nent of the ZDI maps, in which only the most robust large-scale 
information is retained. We find that during the first two epochs 
(2004.6 and 2006.7) the global field is nearly axisymmetric and 
is reminiscent of a dipole aligned with the stellar rotational axis. 
A drastic change of the large-scale magnetic topology occurs be- 
tween epochs 2006.7 and 2007.6. Simultaneously with a sharp 
decline in the total magnetic field energy, the field becomes more 
complex and decidedly non-axisymmetric. 



6. Discussion 

6.1. Comparison with previous ZDI studies 

The only other RSCVn system repeatedly studied with ZDI 
is HR1099 (V71 1Tau). Th e papers by |Donati| ^999\ , |Donati| 
et al.H200"3] ), and |Petit et al.| ( [2004 ) presented magnetic field and 
brightness distributions recovered for about five epochs each, 
spanning the period from 1991 to 2002. In all these studies mag- 
netic field reconstruction was carried out assuming an immac- 
ulate photosphere and the local line profiles were treated with a 
Gaussian approximation or using LSD profiles of slowly rotating 
inactive standards. 

The observed LSD Stokes V profiles of HR 1099 have a typi- 
cal peak-to-peak amplitude of 0.15%, whereas the reconstructed 
field intensities are of the order of a few hundred G on average 
and reach up to ~ 1 kG locally. This is comparable to our obser- 
vational data and inversion results for II Peg. The ZDI studies 
of HR 1099 revealed dominant azimuthal magnetic fields, of- 
ten arranged in unipolar rings encircling the star at a certain 
latitude. Repeated magnetic inversions suggested stability of 
these structures on the time scales of several years. The au- 
thors attributed these horizontal fields to a global toroidal mag- 
netic component produced by a non- solar dynamo mechanism 
distributed throughout the stellar convection zone. 
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Fig. 12. Root-mean- square values of the diff'erent magnetic field com- 
ponents calculated over each ZDI map: a) radial field, b) meridional 
field, c) azimuthal field, d) total field. The relative sizes of symbols and 
their colours correspond to the phase coverage of individual data sets, 
similar to Fig. [3] 



Compared to these studies of HR 1099, our ZDI maps of 
II Peg show a considerably smaller relative contribution of az- 
imuthal fields. We still find a dominant toroidal component; 
however, these results are not fully equivalent nor easily com- 
parable to those of, for example, |Donati| ( |1999 |) because here 
we use a harmonic representation of the magnetic field topology 
and hence are able to disentangle toroidal and poloidal contribu- 
tions to the azimuthal field, whereas previous studies of HR 1099 
completely ignored the poloidal contribution to the azimuthal 
field. 

It is clear that our data contain no evidence of the persis- 
tent unipolar azimuthal ring-like magnetic structures similar to 
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Fig. 13. Long-term evolution of diff'erent contributions to the mag- 
netic field topology of II Peg. The panels show: a) variation of the to- 
tal magnetic energy, b) the relative contributions of the poloidal (filled 
symbols) and toroidal (open symbols) harmonic components, and c) 
the relative contributions of the axisymmetric (filled symbols) and non- 
axisymmetric (open symbols) components. The relative sizes of sym- 
bols and their colours correspond to the phase coverage of individual 
data sets, similar to Fig. [3] 



those reported for HR 1099. Thus, either the dynamo mecha- 
nism operates diflTerently in the two RS CVn stars with nearly 
identical fundamental parameters, or azimuthal fields may rep- 
resent an artefact owing to a simplified ZDI approach adopted 
for HR 1099 ( Carroll et al.||2009al ). To this end, we note that 
any axially- symmetric structure appearing in stellar DI maps 
must be carefully examined to exclude possible systematic bi- 
ases and, at least, must be connected to the stationary features in 
the observed profiles to prove its reality. Neither was done for 
HR 1099. Faced with this discrepancy between major features 
of the magnetic maps of II Peg and HR 1099, we conclude that a 
definite confirmation of the dominant ring-like azimuthal fields 
in the latter star and associated inferences about field genera- 
tion mechanisms must await improved inversion methodologies 
and observations in all four Stokes parameters ( [Kochukhov et al. 
[2QTT] ). 
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Fig. 14. Extended magnetosphere of II Peg determined using potential field extrapolation from the radial component of ZDI maps. The 
star is shown at four equidistant rotation phases (columns) for the four epochs with best phase-coverage (rows). The underlying spherical map 
corresponds to the radial magnetic field component. The open and closed magnetic field lines are shown in diff'erent colour 



6.2. Relation between temperature and magnetic field maps 

Throughout the entire series of ZDI images reconstructed for 
II Peg we see no obvious spatial correlation between the low- 
temperature spots and the strongest magnetic field features. This 
suggests that a significant fraction of magnetic flux is not asso- 
ciated with cool spots. Although a similar conclusion has been 
reached by p revious studies ( [Donati & Collier Cameron||l997 



Donati|1999| ), it was unclear whether this represented a genuine 
characteristic of an active- star stellar surface structure or an in- 
version artefact coming from an inconsistent modelling of the 
magnetic and temperature spots. Our work confirms the lack of 
the field-spot correlation based on a self-consistent and phys- 
ically realistic analysis of the circular p olarisation in spectral 
lines. Using a similar approach, [Carroll et al. ( 2007 ) also failed 
to detect a strong correlation between magnetic and temperature 
features on the surface of II Peg. These results may be inter- 



preted as evidence that current ZDI maps are mostly sensitive 
to magnetic fields at photospheric temperature and entirely miss 
the very strong fields inside cool spots. Instead, numerical ex- 
periments predict that self-consistent ZDI should be capable of 
recovering fields inside cool spots even if the temperature con- 
trast is as large as 1500 K ( [Kochukhov & Piskunov|2009t [Rosen 
|& Kochukhov|2012| ). 

A limited spatial resolution of the ZDI maps may be an- 
other reason for not seeing a link between fields and cool spots. 
Both the magnetic and temperature distributions obtained using 
a spectral inversion technique reveal only the largest- scale struc- 
tures, which are probably not monolithic but consist of many 
smaller spots with diflferent geometries and field polarities. Then 
a local correlation between low-temperature spots and magnetic 
fields may be washed out in the current generation of ZDI maps. 
One can note that the presence of a substantial unresolved small- 
scale magnetic flux implies a much stronger average magnetic 
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field strength than the one inferred from the ZDI analysis of cir- 
cular polarisation. Highly inconsistent results of ZDI and Zee- 
man splitting studies of low-mass stars represent an example of 



this situation for a difiTerent dynamo regime ( [Reiners & Basri 
[2009j . 

At the same time, it is not entirely obvious from a theoreti- 
cal standpoint that a one-to-one relation between the cool spots 
and magnetic fields stemming from the solar paradigm can be 
universally extended to other types of cool active stars. There 
is an increasing amount of theoretical evidence (e.g. |Chan|2QQ7[ 
|Kapyla et al.| 2Ql lb ; M antere et al.|2Ql l| of a pure hydrodynam- 
ical instability leading to the generation of large-scale vortices 
in the rapid rotation regime. These structures have so far been 
found only in local Cartesian simulations of turbulent convec- 
tion, the sizes of the vortices always being very close to the box 
size, suggesting that these structures may have globally signifi- 
cant spatial extents. Depending on the rotation rate, either cool, 
cyclonic vortices for intermediate rotation, or warm anticyclonic 
vortices are excited, the temperature contrast being of the order 
of ten percent. This instability might contribute to the genera- 
tion of magnetic fields independent from temperature structures 
in rapid rotators. So far, however, these structures have been 
detected neither in more realistic spherical geometry nor in the 
magnetohydrodynamic regime (see e.g. |Kapyla et al.|2012b| ). 



6.3. Interpretation in terms of dynamo theory 

The picture arising from an extensive set of previous photomet- 
ric and spectroscopic observations of II Peg suggests that the sur- 
face magnetic field of this object concentrates on one or two ac- 
tive longitudes, i.e. is highly non-axisymmetric, and that these 
active longitudes evolve dynamically over time. During the 
epoch 1994-2002 , a persistent drift of the active longitude has 
been confirmed (Berdyugina et al.|1998a| |1 999[|Lindborg et aL 
20TTI). It is not detectable during 2004-2010 ([Hackman eTaL 



2012| ) which is supported by the ZDI maps presented in this pa- 
per. Furthermore, our new results give an indication that the 
strength of the magnetic field has been monotonically decreas- 
ing, at least during 2004-2009; it seems that a minimum was 
reached at about 2009, after which the magnetic field strength 
started increasing again. The analysis of the energy contained in 
the poloidal and toroidal components as function of time shows 
that while in the beginning of the dataset the field was predom- 
inantly poloidal, the portion of the toroidal field is increasing 
nearly linearly with time, and is dominating at the end of the 
dataset. This is also reflected by the increasing contribution 
of the non-axisymmetric component exceeding the energy con- 
tained in the axisymmetric modes for 2009-2010. All these find- 
ings together hint towards a possible minimum in the star's mag- 
netic activity cycle, during which the magnetic field tends to be 
more poloidal and axisymmetric, accompanied with the signa- 
ture of the non-axisymmetric drifting dynamo wave getting too 
weak to be detectable. Unfortunately, the datasets, especially the 
series of ZDI maps, are too short to make decisive conclusions 
on the cyclic nature of the magnetic field. In any case it seems 
evident that the magnetic field on global scale is far from static. 

How can this be understood in terms of dynamo theory? 
In the solar case, the internal rotation and its non-uniformities 
are known from helioseismic inversions, while in the case of 
other stars, photometric period variations, interpreted as indi- 
rect proxies of stellar surface diflTerential rotation, are normally 
much smaller than the solar value. Theoretically this is conceiv- 
able, as it has been predicted that the faster the star rotates, the 
smaller the non-uniformities in its rotation rate will be (see e.g. 



Kuker & Rudiger|[2005| ). This means that in the rapidly rotat- 
ing late-type stars the operation of the dynamo should rely more 
strongly on the collective inductive action of convective turbu- 
lence (called the a eflTect), and less on the non-uniformities of the 
rotational velocity (called the Q eflTect), the dynamo therefore be- 
ing more of the type than the OfQ solar dynamo. The simplest 
of these systems (see e.g. Krause & Raedler| |1980) excite dy- 
namo modes that are non-axisymmetric but show no oscillations, 
although drifts of the magnetic structure, i.e. azimuthal dynamo 
waves, are typical. As more and more observational evidence 
on dynamically changing magnetic fields is being gathered, it 
has become evident that this simple picture is not adequate. It 
has been suggested that either these objects have more diflTeren- 



tial rotation than predicted by theory (e.g. Elstner & Korhonen 



2005| ), or that the mean-field transport coeflftcients describing 



the convective turbulence are far too simple. Indeed, oscillat- 
ing dynamo solutions in the regime have been found with 
more complex profiles (jB^yshnikova & Shukurov|1987[|Mitra 



|et al.|2010t [Kapyla et al.|2012b| ). We also note that some obser- 
vational evidence exists supporting the idea of stars showing a 
larger amount of diflferential rotation than actually predicted by 
the theoretical models, although the discrepancy between obser- 
vations and theory appears to be quite small ( Hall| 1991 [[Collier 



Cameron 



[2007| ). Direct numerical simulations in spherical ge- 



ometry also show results consistent with the quenching of rel- 
ative diflTerential rotation with increasing rotation rate ( Kapyla| 
et al.|20iri| , although it is still challenging to relate these mod- 



els to real stars. Therefore, we cannot completely rule out the 
existence of enough diflTerential rotation in rapid rotators to be 
significant for the dynamo mechanism. 

For our purposes, it is relevant to compare our observational 
results with direct numerical simulations of turbulent convec- 
tion. Extensive parameter studies have been performed partic- 
ularly in Cartesian geometry ( [Kapyla et al.||2Q12b[ ), while such 
studies in global spherical geometry remain challenging ( Miesch| 



& Toomre^2009 ). One clear shortcoming of the local Cartesian 



models is that the diflTerential rotation cannot self-consistently 
emerge as a result of the modelling but needs to be imposed. The 
global spherical models can grasp this aspect, but only a few are 
successful in reproducing the solar rotation profile (see Mieschl 
et al.||2006 ). The Cartesian studies clearly indicate that oscil- 



latory a^ynamos are quite natural in the rotation-dominated 
regime. Typical solutions (e.g. Kapyla et al.' 2012b) show ra- 
dial and azimuthal components of nearly equal strengths, with a 
rough 7r/2 phase separation in the cycle. Inclusion of shear into 
such a system has two principal efl'ects. Firstly, the azimuthal 
component grows in strength versus the radial component, and 
quite often the components are in anti-phase, the sign changes 
and minima occur simultaneously. The ZDI data indicates that 
the azimuthal field component is not completely synchronised 
with the radial and meridional fields, i.e. they do not seem to 
grow/decline simultaneously. This, again, is more consistent 
with the scenario than the aQ. picture. 

The dominating non-axisymmetric topology, azimuthal dy- 
namo waves, and even the time dependence of the magnetic 
fields in rapid rotators can be quite readily understood in terms of 
dynamo theory relying only on the a mechanism; aQ. dynamos, 
in contrast, produce latitudinal dynamo waves (such as the solar 
butterfly diagram), mostly axisymmetric fields, and oscillatory 
solutions are the preferentially excited ones (see e.g. |Steenbeck[ 
"& Krause 1969 ). The presence of a significant axisymmetric 
contribution, as was found in the case of II Peg in this study, is 
difl&cult to explain with a pure dynamo mechanism. To "ax- 
isymmetrise" some part of the dynamo solution, some differen- 
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tial rotation will probably be needed; therefore our ZDI results 
most reasonably point towards an a^Q. dynamo operating on the 
object under study. 
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Table 1. The journal of spectropolarimetric observations of II Peg. 
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13796 


12.1 


62.5 + 3.7 




09/09/2006 


53987.7055 


0.051 


368 


13427 


13.0 


83.1+3.7 




11/09/2006 


53989.7039 


0.348 


343 


12866 


6.7 


48.7 + 4.0 




12/09/2006 


53990.6496 


0.488 


438 


16330 


5.3 


89.6 + 3.1 




13/09/2006 


53991.6798 


0.642 


185 


7181 


8.2 


4.7 + 7.1 




14/09/2006 


53992.6805 


0.790 


232 


8898 


10.3 


86.6 + 5.7 


2006.9/63% 


01/12/2006 


54071.4660 


0.507 


215 


8461 


8.2 


22.4 + 5.9 




03/12/2006 


54072.5264 


0.665 


162 


5987 


12.2 


-45.0 + 8.2 




04/12/2006 


54074.4128 


0.945 


199 


7634 


14.4 


64.0 + 6.5 




05/12/2006 


54075.4717 


0.103 


172 


6316 


14.3 


-10.8 + 7.8 




06/12/2006 


54076.4425 


0.247 


203 


7614 


11.1 


-63.6 + 6.5 




07/12/2006 


54077.4278 


0.394 


166 


6691 


10.8 


17.3 + 7.5 




08/12/2006 


54078.4148 


0.540 


148 


5822 


9.2 


-65.7 + 8.6 


L\j\j 1 XijoV /C 


1 0/07/0007 


^ZL^OO ^0^1 


O ^07 


ZJ J 


Q7^ZL 
y /D4 


J. 4 


70 7 -u ^ ^ 
— /Z. / ± J.J 




Z\J/\J 1 1 Z\j\j 1 


^ZL^OI 71 O^ 


O 7ZL7 


9Q0 

zyu 


1 OQ^^ 

Wjyjj 


1 9 A 


7ZL -1- ZL 7 
— / +.U ± +. / 




21/07/2007 


54302 6741 


0.891 


299 


11346 


13.8 


33.8 + 4.5 




22/07/2007 


54303.6297 


0.033 


291 


11222 


9.2 


19.5 + 4.6 




23/07/2007 


54304.6726 


0.188 


265 


9659 


10.4 


41.4 + 5.3 




24/07/2007 


54305.7046 


0.341 


325 


12704 


6.7 


59.0 + 4.0 




25/07/2007 


54306.6431 


0.481 


343 


13228 


6.4 


-18.6 + 3.9 




26/07/2007 


54307.7051 


0.639 


266 


10224 


6.5 


-108.0 + 5.0 




27/07/2007 


54308.6486 


0.779 


300 


11357 


15.6 


-81.1+4.5 




28/07/2007 


54309.6977 


0.935 


242 


9075 


10.5 


32.3 + 5.6 


2001.9163% 


22/11/2007 


54427.4366 


0.445 


264 


10021 


7.2 


-47.7 + 5.1 




25/11/2007 


54430.4725 


0.896 


166 


5998 


7.6 


58.1 + 8.1 




26/11/2007 


54431.4755 


0.045 


260 


9451 


16.3 


67.2 + 5.3 




27/11/2007 


54432.4542 


0.191 


378 


10369 


9.2 


59.8 + 5.0 




30/11/2007 


54435.3169 


0.617 


128 


4469 


14.1 


-130.2 + 11.1 




01/12/2007 


54436.4705 


0.788 


326 


V115\ 


7.0 


-21.8 + 4.1 




02/12/2007 


54437.4485 


0.934 


321 


11758 


7.3 


73.5 + 4.3 


2008.7/40% 


08/09/2008 


54717.6059 


0.597 


328 


12065 


8.6 


-45.4 + 4.3 




11/09/2008 


54720.5529 


0.035 


313 


12043 


2.1 


-1.7 + 4.3 




13/09/2008 


54722.6254 


0.343 


348 


12863 


12.0 


34.3 + 4.0 




14/09/2008 


54723.5442 


0.480 


330 


12382 


7.6 


26.5 + 4.2 


2008.9/59% 


08/12/2008 


54809.3926 


0.247 


415 


15398 


7.0 


45.5 + 3.4 
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Table 1. Continued. 



Epoch/ 


UT Date 


HJD 


Rotational 


S/N 


S/N 


LSD y„,ax 


(5z> 


Phase coverage 




-2400000 


phase 


(V) 


(LSD V) 


xlO^ 


(G) 




09/12/2008 


54810.4185 


0.399 


114 


3932 


11.5 


8.6 + 12.2 




10/12/2008 


54811.4780 


0.557 


242 


8764 


10.0 


-41.6 + 5.8 




11/12/2008 


54812.4383 


0.700 


283 


10017 


7.9 


-52.2 + 5.0 




13/12/2008 


54814.4361 


0.997 


368 


13643 


3.1 


34.2 + 3.7 




14/12/2008 


54815.4722 


0.151 


292 


10561 


3.6 


20.1+4.8 


2009.7/87% 


26/08/2009 


55069.6762 


0.954 


328 


12128 


7.0 


1.4 + 4.2 




27/08/2009 


55070.6893 


0.105 


310 


11060 


8.4 


-72.2 + 4.6 




28/08/2009 


55071.7238 


0.259 


200 


7003 


19.8 


-35.4 + 7.2 




29/08/2009 


55072.6965 


0.404 


347 


12219 


13.7 


53.6 + 4.2 




30/08/2009 


55073.7070 


0.554 


202 


7326 


20.1 


67.1+7.0 




01/09/2009 


55075.5550 


0.829 


274 


9898 


13.8 


40.6 + 5.1 




02/09/2009 


55076.6844 


0.997 


393 


15026 


4.3 


-19.5 + 3.4 




03/09/2009 


55077.6378 


0.138 


443 


16263 


10.7 


-49.8 + 3.2 




04/09/2009 


55078.7420 


0.303 


359 


13033 


14.9 


-2.3 + 3.9 




05/09/2009 


55079.6788 


0.442 


372 


13695 


11.2 


36.4 + 3.7 




06/09/2009 


55080.6736 


0.590 


374 


12305 


19.7 


59.9 + 4.2 




07/09/2009 


55081.7424 


0.749 


177 


5654 


11.0 


30.5 + 9.2 


2010.0/68% 


27/12/2009 


55193.3636 


0.348 


296 


10102 


10.9 


3.4 + 5.0 




28/12/2009 


55194.3791 


0.499 


287 


9793 


13.4 


55.5 + 5.1 




29/12/2009 


55195.3684 


0.647 


397 


14431 


7.3 


32.9 + 3.6 




31/12/2009 


55197.3984 


0.948 


397 


14261 


10.8 


-9.4 + 3.6 




01/01/2010 


55198.4116 


0.099 


338 


12572 


5.0 


-57.5 + 4.1 




02/01/2010 


55199.3720 


0.242 


283 


10120 


16.7 


-63.2 + 5.0 




03/01/2010 


55200.3627 


0.389 


230 


8445 


14.4 


-4.6 + 5.9 




04/01/2010 


55201.3844 


0.541 


235 


8599 


13.3 


50.8 + 5.9 



V/Ie (%) 
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Fig. 7. Same as Fig.|5]for epochs 2005.0 and 2005.6. 
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Fig. 8. Same as Fig.|5]for epochs 2005.9 and 2006.9. 
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2007.9 
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Fig. 9. Same as Fig.|5]for epochs 2007.9 and 2008.7. 
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Fig. 10. Same as Fig.[5]for epochs 2008.9 and 2010.0. 



